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Abstract 

(N : 

We analyze the perturbative implications of the most general high derivative approach to 
quantum gravity based on a diffeomorphism invariant local action. In particular, we con- 
sider the super-renormalizable case with a large number of metric derivatives in the action. 
The structure of ultraviolet divergences is analyzed in some detail. We show that they are 
independent of the gauge fixing condition and the choice of field reparametrization. The 
cosmological counterterm is shown to vanish under certain parameter conditions. We elab- 
orate on the unitarity problem of high derivative approaches and the distribution of masses 
of unphysical ghosts. We also discuss the properties of the low energy regime and explore 
the possibility of having a multi-scale gravity with different scaling regimes compatible with 
Einstein gravity at low energies. Finally, we show that the ultraviolet scaling of matter theo- 
ries is not affected by the quantum corrections of high derivative gravity. As a consequence, 
^ . asymptotic freedom is stable under those quantum gravity corrections. 

J3; 

Introduction 

The formulation of a consistent theory of Quantum Gravity is still one of the major 
challenges in theoretical physics. One of the main fundamental problems is that there is no 
experimental evidence of any Quantum Gravity effect Jl|, whereas the classical theory covers 
from cosmology to current precision tests with great success. The situation is considerably 
more dramatic that for Quantum Electrodynamics at the end of the forties. The reason 
being that in the last case, although the theory was not completely consistent, there was 
experimental evidence on the existence of relativistic quantum effects associated to electron 
dynamics. Nature provided confidence on the co-existence of special relativity and quantum 
mechanics although field theory was not yet at sight. An analogous conviction does not exist 
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in Quantum Gravity. On the other hand since the quantum effects can only, in principle, be 
observed at the Plank scale which lies too far from current experiments, the lack of experi- 
mental data is not unnatural. Consequently, any theoretical model of Quantum Gravity has 
to be necessarily based on theoretical principles without phenomenological constraints. 

On spite of this freedom of the theory of Quantum Gravity, the construction of consistent 
models meets very serious difficulties. First, the quantum theory based on the Einstein- 
Hilbert action is non-renormalizable 0, ||. The radiative corrections to the effective action 
contain divergences which involve a number of derivatives of the metric in the counterterms 
which is increasing with the number of loops. Indeed, one can use the general theorems on 
covariant renormalization (see [|J and references therin) to show that all these counterterms 
are general covariant local expressions. Thus, it is possible to remove all the divergences 
starting from the Einstein theory by adding an infinite tower of possible higher derivative 
terms, regarded as the perturbations. In that case the definition of the renormalized theory 
requires, in general, the introduction of an infinite number of renormalization conditions, 
unless some parameter reduction mechanism appears p. Recently, it was pointed out f| that 
some large distances effects remain, however, independent of the renormalization conditions 
involving higher order terms of the effective action. In particular this occurs for the quantum 
corrections to the gravitational Newtonian potential. However, the analysis of other effects 
like the scaling dependence of the cosmological and gravitational constants, depends on the 
concrete formulation of the quantum theory. 

One can always imagine that the extremely high energy UV regime is described by qual- 
itatively different theory like string theory, free of renormalizability problems, which will 
provide a natural reduction parameter scheme at intermediate scales. Having this as per- 
spective, field theoretical approaches to quantum gravity should be considered as effective 
field theories. However the predictions of string theory quantum gravity effects at the in- 
termediate energies still remain unveiled. On the other hand, the string effective action 
is well defined only on shell. Off shell continuations are not uniquely defined because of 
reparametrization invariance [7|— [I5|. For such a reason, we will consider the most general 
effective actionP] which can be generated by string theory [T7 

N+2 . 

Ses(g^) = hm V a 2n m 4 ~ 2n / d A x^g 2n {d\9^v), (1) 

N ^°°n=0 J 

where m 2 is the only dimensional parameter of our fundamental theory (in string theory it 
is the string tension 1/a') while 02n{dxg^ v ) denotes the general covariant scalar terms con- 
taining 2n derivatives of the metric g^ v and the constants a n are dimensionless couplings. 
Although theories with higher derivatives like @) are in general non-unitary at the quan- 
tum level, string theory is both unitary and renormalizable. In particular, one can always 

2 We do not consider the dilaton and the antisymmetric fields for the sake of simplicity throughout this 
paper. The inclusion of such fields does not introduce qualitative changes in the results. 
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choose a special parametrization without unphysical ghosts U , [14| , although from a pure 



string point of view there are no means to distinguish this special parametrization of the 



metric except for the absence of ghosts |J (see also [|16j). However if one considers an ap- 
proximation to the effective theory and makes a truncation of the series ([J), the unitarity 
problem reappears. The problem shows up even at the classical level where, because of the 



existence of the unphysical ghosts with negative energy, it leads to classical instabilities 
The quantum unitarity problem has been considered in great detail in the particular case of 
(four derivative) i? 2 -gravity p|, ||, p|-p|. In that case, the ultraviolet behaviour of the 



propagators and vertices leads to a renormalizable theory |2D, M. Unfortunately the particle 
content of this theory contains, besides the massless graviton, a spin-2 massive unphysical 
ghost which violates unitarity EO, [Tj|. Despite a lot of interesting attempts to solve the 



unitarity problem in i? 2 -gravity |24| it turns out that the massive ghost is not removed 
by radiative corrections, and therefore high derivative theory can not be considered as the 
fundamental theory of quantum gravity f\ The appearance of the massive ghosts is the price 
to pay for renormalizability. Their contributions are essential to reduce the dimension of 
counterterms. However, since the masses of the ghosts are of the order of the Planck mass, 
the fourth derivative quantum gravity can be successfully used as an effective theory at the 
energies below this scale, where the unphysical ghosts are not generated p6|-|f29[ . 



The aim of the present paper is to extend those results for the general case of high 
derivative gravity. We will truncate the effective string theory action S e g taking a finite 
value for N. If N > the corresponding theory contains more that four derivatives of the 
metric and it becomes super-renormalizable. This provides a natural framework to study 
the possibility of different scaling regimes compatible with Einstein gravity at low energies 
and different scenarios for the ultraviolet behaviour of the cosmological constant. 

The organization of the paper is as follows. In section 2 we consider the quantization of 
the general higher derivative theory and show its super-renormalizability. In section 3 the 
general structure of ultraviolet divergences is analyzed and the cosmological counterterm is 
explicitly calculated for the general case. We also study the gauge fixing independence of the 
counterterms. The low energy regime of the theory is examined and shown to be equal to 
that of Einstein gravity in section 4, where we also discuss the spectrum of massive ghosts, 
and the possibilities of having different scaling regimes. In section 5 we study the coupling of 
high derivative gravity to general gauge theories and show that their beta functions are not 
modified by the gravitational corrections. Finally, a discussion of the results is developed in 
section 6, and some technical aspects are postponed to three appendices. 

2. Gauge fixing, quantization and power counting. 

3 It is a remarkable fact that all attempts to derive a Field Theory of Gravity from non-commutative 
geometry lead to i? 2 -gravity 
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The first two terms of the action S e g are the cosmological term 

s §s($iw) = m4 «o / d 4 x^fg 

and the Hilbert-Einstein action, 

The third term (n = 2) leads to the mentioned fourth derivative gravity 

SSW) = j d A x^ (a\ R aP K# + a\ R 2 



(see fl3(J for an introduction and more complete references). Higher order terms (n > 2) can 
be expressed in terms of the Riemann curvature tensor, Ricci tensor, scalar curvature and 
their covariant derivatives. For n = 3 we have two different types of terms R^ and Vi2 • Vi? 
and for n = 4 we have R 4 , RVR-VR, V 2 R-V 2 R. N-th order include terms from R N+2 and 
R N ~ 2 V-R " Vi? " to R,,, □ i? ". The dots indicate all possible contractions of tensor indices. 
Perturbation theory is generated by the standard expansion around the flat metric 

The contributions of the N > terms to the propagator can only come from terms of second 
order in the curvature. Terms of third and higher orders in the curvature contribute only to 
vertices, because they involve more than two fields. Using Bianchi identities 

^xRfiua/3 + V ' pRfj,u\a + Vq-R^/ja = (2) 

one can easily reduce the terms of second order in the curvature (with total 2n derivatives) 
to the form 

a x n R a/3 a n ~ 2 R^ + a 2 n R D n ~ 2 R + a 3 n R, ua/3 n n ~ 2 R^ (3) 
Now one can use (§) again as it was done in || for the N — 1 (six-derivative) terms to get 
R^a n - 2 R^ = -VxR^a^^R^ + 0(R 3 ) + V^ 

= AR ap n n - 2 R a ? -RD n ~ 2 R+0(R 3 ) + (4) 

This is indeed similar to the relation which takes place in the fourth derivative gravity 
because of the topological nature of the Gauss-Bonnet term 0. Of course, for any n ^ 2 
the term 

/ &xyf=g (R^ap n" -2 R^ va(i - 4 R af3 n n " 2 R a/3 + R a n ~ 2 Rj (5) 
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is not topological and gives rise to non-trivial contributions to the vertices. Since the 0(R 3 ) 
terms do not contribute to the propagator, the relevant terms of order n > 2 can be written 

as 

a\ R af3 U n ~ 2 R a ? + a 2 n R U n ~ 2 R (6) 

If ot\ N+i 7^ and ct\ N+A + 30^+4 ^ 0, the highest order terms of (|6]) are nondegenerate 
once we introduce a gauge fixing and asymptotically they behave like 0{k 2n+i ) for large 
momenta. The gauge fixing condition can be introduced within the standard Faddeev- 
Popov prescription. In order to improve the regularity properties of the quantum fields it is 
convenient to add a higher derivative covariant operator for the longitudinal modes in the 
gauge fixing condition [3(J (see also |H|]-|33|] for a general discussion on gauge theories) 



where 



Sgf = / d^x X , C» X u (7) 



X» = d x hl-(3d»h x 



A 



N+2 



□ 



^ = "E kn^n + ^n-l^VH — (8) 

n=2 X - m ' 

P,<^n,Jn being dimensionless gauge fixing parameters. Regardless to the number of deriva- 
tives, the propagator of the quantum metric can be written in terms of the irreducible spin 
2,1,0 projectors and some spin- zero transfer operators pCjj . By choosing the gauge fixing 
parameters (3,a n ,^ n in a special way one can always remove all the spin-1 states from the 
spectrum. The spin-2 states are gauge- fixing independent, just as in the fourth derivative 
gravity, and if a\ N+A ^ 0, the propagator of the spin two states has the asymptotic ultravi- 
olet behaviour 0{k^ 2N+ ^). The same occurs for the propagator of spin-0 states but there 
the condition reads ct\ N+i + 3a| Ar+4 ^ 0. Throughout we shall assume that both conditions 
are satisfied and therefore the propagator of the quantum metric behaves like 0(fc~( 2Ar+4 < ) ) 
for large momenta. One can absorb the determinant of the covariance operator C^ v into the 



action of Faddeev- Popov ghosts ||2q| , 



>FP 



J d A x^ [ c x V u ( V^c A + Va^) - 2 p c x C Xu V u (V„c")] (9) 



and get that their propagator also behaves like 0(k~( 2N+4 ^) in the UV regime. The extra 
contribution of the operator C into the ghost sector has to be compensated by the corre- 
sponding power of the determinant detC Ml/ |2(| (see also [3lH— 1|33|| for a general discussion 
on gauge theories). In summary, the partition function is given by 

[<*V] [6Cr] [5c X ] (detCT* jSM+iSg+iSFr ( 1Q ) 
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With that gauge fixing choice the propagators of both metric and ghosts have the same 
asymptotic behaviour in the UV limit 0{k~( 2N+/ ^). Unfortunately the interaction vertices 
pick up all possible number of derivatives from zero to 2N + 4, in such a way that one loop 
divergences remain unregularized. Indeed, if we evaluate the superficial degree of divergency 
D of an arbitrary p-loop diagram with ri2 r vertices with 2r derivatives r — 0, 1, • • - , N + 2 
we get 

AT+l 

D + d ext =A+2N-2Np- ^(2A^ + 4-2r) n 2r (11) 

r=0 

where d ext is the total number of derivatives acting on external lines. From the power 
counting identity (|TT| ) it follows that divergent graphs satisfy the inequality 

d ext <A + 2N-2Np (12) 

which implies that divergences can only appear for higher loops diagrams with iV = 0, 1, 2. 
For N > 2 divergences appear only in the one-loop diagrams. For any N > the divergent 
terms involve less powers of the curvature than those bare action, and moreover, this power 
is decreasing with the loop order. We remark that the covariance of all the counterterms (in 
an invariant regularization) is guaranteed by the general theorems Q which can be trivially 
generalized for diffeomorphism invariant theories of the form S e g. 

Thus, for any finite N > 0, the structure of divergences becomes simpler, the theory is 
super-renormalizable. For N > 2 the situation is even better because all the divergences 
appear at one loop order of perturbation theory. Now, in any case the structure of one 
loop divergences is not changed. This is very similar to what happens in gauge theories 
where the method of high covariant derivatives does not smooths the behaviour of one loop 



contributions pifl . 



If iV > 2, since the only UV divergences appear in one loop diagrams, they can be removed 
by Pauli-Villars determinant regulators using the methods first introduced for gauge theories. 
In such a case we end up with a completely finite theory (all beta functions vanish) and the 
corresponding theory may be considered as a complete regularization of Quantum Gravity. 
This is nothing but the implementation for quantum gravity of the high derivative method 
introduced by Slavnov |[34!| . In such a regularization one can look for non-perturbative 
effects and specially for non-trivial fixed points where the theory might be unitary and 
renormalizable. If such fixed points would exist most of the technical problems mentioned 
in the introduction will be overcome. 

For any N > the theory is super-renormalizable, because the local covariant countert- 
erms have less derivatives than the classical action and the coefficients of the terms with 
more derivatives do not need any kind of infinite renormalization. Thus, the bare values of 
those coefficients can be kept finite. Their explicit value can only be fixed by looking at the 
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phenomenological implications of the theory or the predictions of the fundamental theory 
which generates our theory as low energy effective theory. The only sector of the theory 
which is subject of infinite renormalization is of the form 

J d^x-s/^g | a\ R a /3 R af3 + a\ R 2 + «2 m 2 R + a m 4 | + (surface terms) (13) 

After renormalization the scaling behaviour of the physical quantities is governed by a renor- 
malization group equation. Since the only parameters which undergo an infinite renormal- 
ization are those contained in the lower derivative part (|T3|) of the action, the scale depen- 
dence of the effective theory S e s is encoded in the beta-functions for the four parameters 
j3 a i, /3 a 2, p a2 , /3 ao . We shall focus on the evaluation of these beta- functions. 

Indeed for large N the Feynman rules are rather involved and in general even one loop 
calculations are very difficult. In the next section we develop a method to calculate these 
beta-functions in one-loop approximation for the most general case, and we shall perform 
an explicit calculation of the cosmological constant beta function (3^ = f3 ao . Our method is 
essentially based on the techniques introduced by Barvinsky and Vilkovisky ||37|| , which can 
be, in principle, applied for the calculation of (3 a i, (3 a 2, f3 a2 as well. The only difficulty is 
that the manipulation of algebraic expressions becomes much more involved in those cases. 



3. One-loop results. 

Before entering into the calculation of beta-functions, let us show that they are indepen- 
dent of the choice of the gauge fixing parameters (|8|) and the parametrization of quantum 
gravitational fields. Local counterterms can be, in principle, gauge and parametrization de- 
pendent, but we will show that these dependence vanishes for high derivative gravity. This 
follows from the fact that for N > 2 the only UV divergences appear at one loop order and 
the explicit relation which exists between the divergent one-loop counter-terms associated to 
two different sets of the gauge fixing parameters o = (cr n ,P,j n ) and a = (a n ,^,%) |28|, [35[ 
(see also Appendix B for a selfcontained derivation). 

r 5S 

T d iv(o r ) - r div (a) = / d 4 xy^g - S^(^ a/S , a, a), (14) 

^fiuido-p, o", o") being some local function of metric and gauge fixing parameters. Power count- 
ing tells us that for any choice of the gauge fixing parameters the divergent counterterms 
are local expressions with up to four derivatives of the metric. Thus, the left hand side of 
the identity (|i4|) can only contain such a type of terms. But the right hand side which is 
proportional to classical motion equations contains terms with 2N + 4 derivatives. Hence the 
equality can hold if only and only if T,^ u (g a/3 ,a,a) is identically zero. One-loop divergences 
are, therefore, gauge independent. Since N > 2 there are not more divergences in higher 
loops and, therefore, we have proved that all the beta- functions /3 a i, f3 a 2, (3 a2 , [3 ao do not 



7 



depend on the gauge fixing condition. In fact, by means of Nielsen identities [g(J it can be 
proved that even in the case < N < 2 there is no dependence on the gauge fixing param- 
eters in the divergent counter-terms generated by higher loops. The independence on the 
parametrization of the quantum metric can be proved in a similar way because two different 
reparametrizations lead to divergent contributions which satisfy an equation similar to ( ^4|) . 

Let us now calculate the one-loop radiative corrections. The leading ultraviolet term is 
given by 

a\ N+A R al3 D N R a ? + a 2 2N+4 R D N R (15) 

We shall assume that ot\ N+i ^ and a\ N+A + 3o3, N+A ^ 0. In the background field 
method the metric is split into a background metric g^ v and quantum metric h^, 

We introduce the background field gauge fixing condition 

Sgf = Jd'x^9 X, C" 1 " Xu (16) 

where 

C^ = -- (^□ + 7 V' i V 1 '- VV) (-^Y (17) 
a \m 2 ) 

which is a covariant generalization of (0) defined by replacing ordinary derivatives by co- 
variant derivatives with respect to the background metric g^ v . The coefficients of the gauge 
fixing condition a, [3, 7 are arbitrary parameters. Many of the coefficients of (|j) have been 
set equal to zero to simplify the calculations, but the divergences will not depend on their 
values. The one-loop effective action is 

r« = % - Tr In H mour - % Tr In M° a - % - Tr In (18) 



where 



5h aa 5h^ v 



+ ^XA_ C \r S Xr 



(19) 

h=0 



and M° is the operator of the FP ghost action 

= + V a V CT -2/3 V CT V a . (20) 
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The contributions of the operators M° and C^ v to the divergences of the effective action 
(|18|) can be easily evaluated by means of Barvinsky-Vilkovisky techniques [0. The result is 



i Tr In M" 




R T \R 




(21) 



and, 



-Tr \nC u u 
2 



iJV 



— Tr In □) - - Tr In (g^ □ + 7 V M V 1/ - W v V A 



1 

Ye 



d Xi 



AN 7 
'~5~ + 10 



R T \R rX + 



IN 



( 20 



3 

20 



R 2 



(22) 



We remark that the last expression does not depend on 7 due to the cancellation first pointed 
out in Ref. |26| . 

The main technical problem is the calculation of Tr InH^ ^. if is a differential operator 
of (2N + 4) order, with coefficients depending on the curvature tensor of the background 
metric and its derivatives, the parameters of the gravitational action, and the gauge fixing 
parameters a, (3, 7. In spite of these difficulties one can derive a general formula for the 
divergent part of Tr In H^^p and then perform an explicit calculation for the cosmological 
counterterm for arbitrary value N. Let us introduce dimensionful coupling constants 



to 



m 



-2n 



a 



2n+4 



by absorbing the mass parameter in the dimensionless couplings a l n of S e s. Since the di- 
vergences do not depend on the explicit values of the gauge fixing parameters, it is very 
convenient to fix their values to simplify the calculations. A convenient choice is 



a 



2u 2 N 



LO 



1 ' 



7 = -^T, P- 
For this choice of the gauge fixing, the operator becomes 



LO 



N 



Au 2 N 



+ 1. 



(23) 



H, 



LO 



piv,af3 



I; p° 
4 \ w % 



LO 



Loir + AlO 



Nj 



16lo 2 n 



9^ 9 



pa 



^\$pa,a(: 



□ 



N+2 



1 1/ AiA2---A2jv+2 V7 V7 
~r V po,a0 V Ai V 



Ai V A 2 • * • 



7, _|_ W TlT 2 -T 2N+ iy y 

A2iV+2 ' vv pa,a0 v T1 v T 



...v 



T~2]V+1 



+u, 



pa,af3 * v\ * 1 



■ v^ + ^v 2 ^- 1 )} 



(24) 



where V, W, U depend on the dimensionful ratios like oo^^/lo^ and on the curvature tensor 
of the background metric and its covariant derivatives. Now we can use the Barvinsky- 
Vilkovisky method f37j to extract the divergences of Tr hiH^^p. The pre-factor of 
does not contribute to the divergences and therefore it can be omitted. Hence 

Tr ]nH maf) = (N + 2) Tr In [8^ aP □] + Tr V^ iajB AlAa "' Aaw + a V Al V Aa • • • V AaiV+2 
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+ TrW^ a ^- T ^V Tl V T2 ■ ■ ■ V T2JV+1 ^ + TrU^ a ^-^V Vl V V2 ■ ■ ■ V„ 2JV ^ 



9 M^.po- » ai » a 2 jv +2 |— |7v+2 0< Pi"P A i A 2JV+2 qJV+2 

+ (terms of higher background dimension) (25) 

Indeed by dimensional arguments the last higher background dimension terms do not con- 
tribute to the divergences. The only term which can not be directly handled by the method 
of p?fl, is that which involves two V matrices. In particular there are new terms coming 



from the commutation of V with the covariant derivatives and the operators d -1 . However, 
since V has dimensions of curvature, and derivatives of V increases its dimension, those 
terms do not contribute to the divergent part of the effective action. Therefore as far as one 
loop divergences are concerned V may be considered as a constant and we can ignore its 
derivatives. The corresponding contribution can be written as Q 

" \ V* M ^ W Ar " W W 1 -^ 2 V* ■ • • V w V*i • • • Vy 2N+2 ^- 4 (26) 

Now all the terms in the expansion admit the direct substitution of the universal trace formu- 
lae of Barvinsky-Vilkovisky and we can derive a general formula for the one-loop divergences 
of the theory. To apply this formula for a general N one needs to expand the action up to 
second order in quantum fields, keeping the first and second orders in background curvature 
while neglecting the higher orders and derivatives of the curvature. The algebra involved 
in the calculation is quite cumbersome. We will restrict ourselves to the calculation of the 
divergent contributions to the cosmological constant counterterm. In order to simplify the 
calculations we remark that the background field method gives the result which is valid for 
any background metric, including the flat one. Thus the cosmological counterterm can be 
derived for the flat background metric, which considerably simplifies the calculation. Ac- 
cording to (|18|) and (|l9| ) the one-loop divergences are related with the Hessian of the classical 



action. If the background metric is flat, one can, therefore, ignore all the terms which have 
more than two powers of curvature, and calculate the divergences of the theory 



r r N 

S=ld 4 x^/=g~ □"fl"" + ■ (27) 

Ln=0 



The cosmological constant term u^^—g does not give divergent contribution for any N > 0, 
and Einstein term is relevant for N — 1 only. Let us first consider the case of iV > 1. The 
relevant part of operator H for the divergent contribution has the form 

iW/? = ^ aN+2 + V^ aP XTpu n N - 1 v x V T W p W a + ^, a/3 Ar ^n iV - 2 V A V T V p V (7 (28) 



4 By the same dimensional reasons the matrix W of (|24| ) can not generate divergences and for such a 
reason we shall not consider its contribution. 
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Using well known expansions of Rn U we obtain the following expressions for V 



U) N LO N \LO N -+- oUJ N ) 

2(cj^_ 1 + 2uo 2 N _ l ) A u) 2 n u n _ x 
H i o ma po 1 H —( — I 



,1 , ,2 



iO 

N-l „At 



U N U N-l r Ar per 

I / I i o 2 N gpv<->a.ft y 

vh(uh + 3u z N ) 



2uo 



N 



Uff-1 + AuJ N-l 



iO 



9\r9a(i^v P " 



N 



J up 



(29) 



and for U 



[Upu,a/3 

4 



Arpo- _ UJ N-2 
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— 0fiv,ap9 9 H —< I " 



N 



2 N 9nu9af39 9" 



2(co N _ 2 + 2u) 2 N _ 2 



to 



N 



lo l n (lo l n + 3uo N 

2 1 12 
r X^,pa , UJ N UJ N~2 ~ ^N^N-2 „ x At „ 

Ow,apO H i 7 1 - o, ,2 \ <W. 9 

N\N * M NJ 



per 



CO 



N-2 At 



2u 

1 gxrgapOvp. 



N 



pa 



UJ 



N 



(30) 



Using the Barvinsky-Vilkovisky trace formulae in Q25D, for flat background metric, we get 
the following expression for the divergences 



Tr \nH, 



pv,af3 



2 : 

\2 



\ J d^xV^g tr (-1 U^gZa + ^ V»"°VW gS P a,rya) , (31) 



where e = (47r) (n — 4) is the parameter of the dimensional regularization, and |37 



( 2 ) i 

g\T P a — 9\t ga + g\p g T a + g\a g P r 



(4) 

g\T\>T>pap>a> = fi'AT gx'r' g P a g P >a> + permutations. 



Now, taking into account (fffi)- 
flat background. The result is 



we calculate the one-loop divergences for the theory on 



oojf(ujlf + 310%) 



u(N) 



K) 2 K + 3^) 2 



v(N) 



(32) 



where 



u(N, N > 1) = (6ulf uj l N _ 2 + lhu? N oj 1 n _ 2 + 3u 2 N _ 2 u^) 
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v(N) = K) 2 K_ 1 + 3^_ 1 ) 2 + 5K_ 1 ) 2 K + 3^) 2 . 



(33) 



For the case N — 1 one needs to change the expression for U in (pO] 



U, 



Xrpa 



(N 



At 1 r ^ Wj + 2wf 



2^(^ + 3^) 9tu,6a ^ u\ 9ap6 ' 



pa 
vp 



Then, after inserting this expression into the formula 
different coefficient u(N = 1). 



(34) 



11) we also obtain (|32!) but with a 



u(N 



2 



3 u\ + 10 uj\ 



(35) 



This formulae ( p2[) (|33| ) and (0) complete the calculation of the cosmological constant 
counterterm in an effective gravity theory S e g. For N > 2 the above expressions are exact 
since there are no any additional divergences for the cosmological constant at higher loops. 
It is important to recall that the result is independent on the choice of the gauge fixing 
condition and on the parametrization of the quantum field. In this respect, it differs from the 



counterterms which appear in quantum gravity with two derivatives |38| and four derivatives 
p6| , f2~8fl . Thus adding the higher derivative terms to the effective theory introduces some 



relevant difference from this viewpoint. 

It should be also interesting to calculate the divergence of the coefficient of the Einstein 
term R, because the corresponding beta function /3q = l3u_i describes the running of the 
gravitational constant with the change of energy scale. The formulas ( |25|) and (|26| ) enable 
one to perform such a calculation for any finite N, but it would require an effort beyond 
our present possibilities. By dimensional arguments the counterterm which is linear in 

CjV-i- We know that all 



curvature has to be a linear combination of the coefficients u>]^_ 1 , 



the coefficients entering in this linear combination will give a factor uj]^ (oj^ + 3u;^-) _:L , but 
the exact value requires explicit calculation. The renormalization of the parameters Wq,0Jq 
of the four derivative terms fll3|) can be also obtained from ( pop and ( |26"D . They are bilinear 



in 



and contain the factor (u> 



uj 1 n + 3lu 



The logarithmic divergent contribution to the renormalization of the cosmological con- 
stant (|32|) yields a renormalization group equation which governs its dependence on the 
energy scale. The corresponding beta-function is 



1 



(4vr) J 



2u(N) 



+ 



v{N) 



;^) 2 K+3^) ; 



(36) 
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One interesting consequence of this general analysis is that it provides different scaling 
scenarios for the running of the cosmological constant (as well as for the other relevant 
couplings). In principle it is possible to impose more constraints on those coefficients to 
suppress all logarithmically divergent contributions. In such a case the corresponding theory 
is finite and all beta functions vanish. But this behaviour only holds for energies higher 
than the effective scale of the theory. When the energy moves to lower values the leading 
ultraviolet terms with 2N + 4 derivatives become irrelevant and new scaling scenarios fl36|) 
governed by the terms with 2N+2 derivatives emerge. This can be seen from our calculations 
taking the limit when the leading ultraviolet coefficients tend to zero. We recover in such 
a case the beta function corresponding to the subleading terms. In this way the theory 
provides a hierarchy of different scaling scenarios from very short distances to cosmological 
distances governed by different beta functions. However, properly speaking to implement 
such multiple scale scenario we will need a set of different dimensionful parameter scales 
Hi, ■ ■ ■ , fj,N. One interesting feature of this scenario is the existence of crossover between 
the different scaling windows appears. The renormalization group flow runs from ultraviolet 
fixed points to infrared fixed points step by step in a continuous non-linear way. 

In this sense an open approach to quantum gravity it is possible, although the com- 
patibility with closed approaches coming from some fundamental theories is also possible by 
imposing the corresponding constraints for the different parameter of the effective action S e g. 
In particular, if a non-trivial fixed point is found the possibility of having a non-perturbative 
approach to quantum gravity is open. The theory defined by scaling limit around this point 
could have unexpected non-perturbative effects. 

4. Structure of mass poles 

In the previous section we have seen that the high derivative theory provides an interesting 
smooth approach to quantum gravity. In this framework the effective theory can exhibit 
interesting features like the vanishing of the coupling constant in the ultraviolet regime. 
Moreover, the theory provides a series of different scaling scenarios for all relevant physical 
quantities governed by different renormalization group fixed points. 

However, as it is well known such a versatility of high derivatives theories implies serious 
problems from an unitarity viewpoint. In particular, the theory contains a plethora of 
unphysical massive ghosts. The analysis of this problem in the framework of string effective 
models |7], § have shown that the massive ghosts disappear from the spectrum in a special 
parametrization of the metric g^ v (the background metric in its target space), while from 
string theory point of view this parametrization does not differ from the others || . 

Recently, it was suggested that in theories with more than two derivatives of the metric 
it could be possible to find a ghost scenario with only one ghost which would be the most 
massive fundamental particle of the theory having negative norm states in the Hilbert space, 
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whereas the other lighter fundamental particles would have positive norm states which are 
compatible with unitarity fl39| . In fact, in the theories of the type considered in the present 
paper the structure of the (euclidean) propagator in the spin-2 and spin-0 reads 

G(Jfe) = (W+4 k 2N+i + l 2N+2 k 2N+2 + l 2N k 2N + --- + l 2 k 2 ) ^ (37) 

where U are real numbers related with the values of uj, uf coefficients in (p7| ). The expression 
fl3?D can be decomposed in terms of simple propagators as 

C(h\ — ^° i , ^2 A N+1 

1 } ~ k 2 k 2 + m 2 + k 2 + m 2+ " " + k 2 + m 2 N+1 1 ' 

where the masses m 2 can be real or complex depending on the values of the coefficients 1% 
The complex masses are always grouped in conjugate pairs. The idea introduced in 



Ref. 3S was based on the assumption that for some physical values of U all the mass poles 



mj are real and positive, < mf < m 2 , < < ■ ■ ■ < m N+l , and the coefficients A r are all 
positive A r > for r = 0, 1, • • • , N, except the last one A N+ i which is negative A N+ i < 0. 
In this case only the heaviest particle would be an unphysical ghost, whereas all the others 
would be ordinary massive particles with positive energy. Thus, the effective theory would 
be unitary only till energies of the order of the mass of the heaviest unphysical particle m 2 N+l , 
although this mass can be chosen much larger than the masses of other particles. Beyond 
such a scale the theory becomes unphysical. Unitarity can only be completely restored for 
all energy scales in the ultimate fundamental theory. 

Unfortunately, this scheme can not work, although from a theoretical viewpoint was very 
appealing. The problem is that the assumptions concerning the behaviour of the massive 
poles m 2 , and their residua Aj in the high derivative theory are not consistent. The reason 
is that (see Appendix A for details) for any real monotone sequence of masses < m\ < 
m 2 < m 3 < • • • < m jv+D the signs of the corresponding residua alternate, i.e. sign [Aj] = 
— sign [A7+1] . This can be easily understood from the the basic property of any polynomial 
with real coefficients and real zeros which establishes that the signs of the slopes at the zeros 
do always alternate between consecutive zeros. Hence the physical assumption made in Ref. 
p9| is never satisfied for real masses. In the case of complex poles we have an even more 
pathological situation. Complex masses m 2 < lead to pairs of unphysical tachyons and 
for the remaining real masses the above argument show that they correspond to alternating 
pairs of particles and ghosts. However, this pathological ghost masses distribution does not 
exclude the existence of a spectrum changing field reparametrization mapping the theory 



into an unitary theory 40 



In any case the dynamical role of unphysical massive ghosts only appears at energies of 
Planck order [p0| . At very low energies (corresponding to the macroscopic length) they do 
not propagate, and the only relevant excitation is that related with the massless graviton. 
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In particular, all infrared quantum effects of the high derivative theory are the same as in 
Einstein gravity. This is specially relevant for quantities with singular infrared behaviour, 
e.g. long range correlation functions. For instance, if we evaluate the one loop quantum 
corrections to the gravitational potential, we get the same result as Donoghue 0, indepen- 
dently of the details of the higher derivative theory we consider. This remarkable result is 
due to the universality of the non-local singular infrared contributions to correlation func- 
tions [41j . Now, in the pure Einstein theory due to the absence of a mass term this singular 
IR behaviour extends to the ultraviolet and in fact can be easily computed in the UV regime, 
where it is linked with the divergent behaviour of the counterterms. The relation is similar 
to the one which exists between the UV counterterms and the beta functions of the coupling 
constant. Now, the UV behaviour of the higher derivative gravity is completely different 
and it is obviously dependent on the details (coefficients and number of derivatives) of the 
effective action. The renormalization group has in such a case a non-linear behaviour which 
goes from the particular UV behaviour associated to our regularized theory to the universal 
infrared Einstein regime. Between these two regimes we have momenta domains where new 
regimes with different scaling behaviour appear provided the masses of the different massive 
ghosts of the theory are widely separated. The nonlinear behaviour of the renormalization 
group allows to go from the UV to the infrared through this series of unstable regimes. The 
relevant result is that the long distance behaviour of the theory is the same as Einstein 
gravity with some small quantum corrections 0. Indeed the contributions of the higher 
derivative terms become more and more relevant as we increase the range of energies or 
what it is the same we go to shorter distances. One of their effects is the generation of a 
running of the gravitational constant with the change of energy scale, and a running of the 
cosmological constant which can become stable at short distances due to the vanishing of 
the corresponding counterterms. 



5. Interaction with matter fields 

Let us now explore the interaction of high derivative gravity with matter and gauge fields. 
In particular, it is interesting to analyze how the quantum effects of this theory can change 
the scaling properties of the matter field correlation functions and the effective potential of 
the Higgs fields. For that purpose we introduce a general gauge model containing spinor, 
vector and scalar Higgs fields with gauge, Yukawa and scalar interactions. The total action 
has the form 

•Stotal = S'cfT + Sgauge (39) 

where S e ff denotes the action of our high derivative gravity (|Ij) and S gaugc is the action of 
the gauge model. We assume, as in the previous sections, that the gravitational propagator 
has an ultraviolet behaviour of order O [ k~( 2N+ ^). 
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The diagrams contributing to any correlation function with external matter lines split 
into two sets. Diagrams with only matter internal lines, and diagrams with an internal 
gravitational propagators. The contributions of the diagrams of the first type (e.g. diagrams 
of the Figure 1) are equal to the contributions of the gauge model in a classical gravitational 
background. It can be shown [|30| (and references therein) that in general the corresponding 
beta-functions for the matter fields couplings (gauge, Yukawa and scalar) are exactly the 
same as in flat space-time. The novel aspects are the need of the non-minimal coupling ^Rcf 1 
for every scalar field 0, the renormalization of £ and the appearance of the corresponding 
beta-function (3$ p0| . 

Indeed, there also appear additional contributions to the renormalization of the param- 
eters c<Jq, lUq, a>_i, cu_2 generated by the loops of matter fields in the diagrams with only 
gravitational external lines (e.g. diagram (2) of Figure 1). Diagrams involving only internal 
and external gravitional lines are only divergent when the inequality ([I2|) holds. Both types 
of divergent counter-terms are of the form (13). 



Let us now consider the most general diagrams which have some external matter lines 
and any kind of internal lines (Figure 2). Here one has to consider the cases N > 1 and 
N = 1 separately. For N > 1 all such diagrams are finite by power counting, thus they do 
not contribute to the beta-functions for the matter field couplings (gauge, Yukawa, scalar), 
which are actually the same as in the theory without quantum gravity. For gauge coupling 
constants the result is even more general. Their beta-function are not affected by any kind 
of the gravitational interaction. In the Einstein gravity it was shown in Ref. ||. The effect 
is essentially based on the non-renormalizability of the Einstein-Maxwell or Einstein- Yang- 
Mills system (power counting does not permit any gravitational contribution to {F® y ) 2 ). In 
the fourth derivative quantum gravity the gravity-independence of the gauge coupling has 



been found in Ref. [|26| as a result of unexpected cancellation of two diagrams. What we 
have shown is that the same result holds for the general theory ([!]) with N > 0. In summary, 
asymptotic freedom is stable under quantization of gravity. 

In the case N = 1, there appear two new divergent diagrams with external scalar lines 
and internal graviton loops (see Figure 3). The divergence of these diagrams is generated by 
the nonminimal interaction term in the scalar (Higgs) sector of the gauge model. In fact the 
divergences generated by such mixed diagrams are always universal, i.e. they do not depend 
on the gauge group or multiplet composition. Consider, for simplicity, one real scalar field. 

1 C 7> ,.,2 / ' 



S sc = Jd^x^ ^_^^0^0 + -m i i + -£ J R0 i -^0 4 J (40) 

Diagrams of Figure 3 can generate UV divergences only if the two derivatives introduced by 
the scalar curvature in | £ R <p 2 term act to the internal gravitational lines. Therefore they 
can only affect the renormalization of the mass m 2 but none of the coupling constants. 
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The corresponding mass counterterm is (see Appendix C) 



This divergence generates a gravitational correction to the known beta-function for the scalar 
mass, which appears due to the gauge, Yukawa and scalar interactions 



,qtotal aknown , rgrav n 

Pm 2 — Pm 2 + Pm' 



8^3 , 1 3 c (c 

Pm 16 7T 2 {u>l + 3u$)* y uj\ 

In summary, in the framework of the higher derivative model of quantum gravity the 
beta functions for the gauge, Yukawa and scalar self-coupling are not affected by quantum 
gravity corrections. The only quantum effect is the correction to the mass beta-function of 
the scalar field (|4lD which only appears for the case N = 1. 

We remark that in fourth derivative quantum gravity the contributions of gravitational 
loops to the beta-functions for the Yukawa and scalar couplings are nontrivial P9l , they can 
even change the asymptotic behaviour of such theories (see also |[30| ). In contrast, for N > 
there is no divergent contribution of gravity to the matter or gauge fields sector. 




6. Conclusions 

The above results show that in the high derivative approach to quantum gravity the UV 
behaviour is so smooth that its corrections to gauge theories do not modify the scaling UV 
behaviour of gauge and matter correlation functions for N > 0. However, in general, the 
pure gravitational selfinteraction is not completely self-regularized. As we have shown one- 
loop divergences are generically of the same type (p~3|) that those of the Einstein theory. This 
is a general feature of non-abelian gauge theories. The analysis of this one-loop divergences 
shows that they can be cancelled by any of the following prescriptions: 

• i) Pauli-Villars covariant regularization 

• ii) Fine tuning of the high derivative couplings 

The theory with a finite number of couplings can be considered either as a regularization 
of quantum gravity or better as an effective theory of quantum gravity. As effective theory 
some of its infrared properties, like the corrections to the Newtonian potential, are universal 
and do not depend on the UV behaviour of the theory. We have shown that the structure of 
mass poles is such that the ghost masses are intercalated between those of particles in such 
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a way that does not allow to recover unitarity unless we restrict to energies below the first 
mass scale of the theory. 

Another interesting feature of the high derivative approach is that the UV scaling regime 
is independent of the gauge parameter and the parametrization of metric fields. 

In the cases where the theory is finite for all momentum scales it might be possible, in 
principle, to find interesting non-perturbative scaling regimes which shall give rise to novel 
non-perturbative approaches to quantum gravity. 

In summary, the high derivative approach offers a good selfconsistent framework to deal 
with quantum gravity effects for intermediate scales between low energy phenomenology and 
high energy fundamental unified theories. 
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Appendix A 

Let us prove now the intercalation property of particles an ghost masses which have been 
used in the main text. We consider for simplicity the euclidean formalism. The euclidean 
propagator G{k 2 ) has the form 



G(k 2 



r 1 



(Al) 



where the masses rrij, j = • ■ ■ N + 1 are functions of the coefficients hj+2] and can be 
rewritten in the form (|38|). Therefore, G(k) has simple poles in k 2 = —m 2 , j — • • • N + 1 
if all the masses are different. 

In the complex plane C, G(z) is a meromorphic function with A^ + 2 simple poles. Let Tj 
be a closed path in C around the pole — m 2 not encircling any other pole. Then, from (A.l), 
we have, for j = ■ ■ ■ N + 1 



I G(z) 



27rzRes 



r 1 

l 2AT+4 



-ITIa 



2-Kil 



-i 



N+l 

2 A' -hi 1 I 

i=0,i^j 



(A2) 



and from 



G(z) = 2mRes 



■N+l 

E 



o (z + m 1Y 3 



2iriAj , 



j = ■ ■ ■ N + 1 



(A3) 
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Then, from (A. 2) and (A. 3) we get 



N+l i 

A* = t£-H II (m 2_ m 2 V j = 0---N + l (A4) 

i=0,i^j \" L i " L j> 

When the masses rrij are real or purely imaginary, their squares m| are real, and once they 
are ordered (m 2 < mf +l Vi = • • • N), (A. 4) tells us that the signs of the residua alternate, 
i.e. aga[Aj] = (-l)^sign[^ +4 ], j — ■ ■ ■ N + 1. 

This result can also be understood in pure algebraic terms. If the polynomial G{k 2 )^ 1 of 
k 2 which defines the propagator only has non-degenerate real zeros (and this is the assump- 
tion we have made) the signs of the slopes at two consecutive zeros alternates. 

On the other hand, if one m 2 is complex, there is another one mf which is its complex 
conjugate, mf = m|*. We can easily show that their residua Aj, Ai are also pairs of complex 
conjugate numbers Aj = A* , whereas the residua Aj corresponding to the remaining real 
masses m 2 have alternating signs. 

Appendix B 

Let us prove in detail the gauge fixing independence of the one-loop divergences by using 
the method of |J7| (see also [|28|, |35| for the fourth derivative gravity). 

Diffeomorphism invariance implies that the classical action is invariant under the in- 
finitesimal transformation 

Sg^ = 2V,a£ a , 2V,a = - (gn\Vv + ^aV m ) 
which leads to the Noether identity 

The gauge fixing term (|7|) introduces a number of gauge parameters <7j,7i,/3. Let us gener- 
ically denote them by a. The propagators G fJjU ' al3 , fl a p of the metric field h^ u and of the 
Faddeev-Popov ghosts c$ x , c^ 2 are defined as usual by the relations 

iW G ^ PU = V 7 , M % & x = 6% (B.2) 

where 

S 2 S Sx P rpa 5xa 



and 



og^v 
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From (B.l) and (B.3) it follows the Ward identity which links the two propagators and 
the gauge fixing functional 

r> o> 3 _ r — r dV\ Tr/ 5S (e>a\ 

og pa og pa 5g X T 

The one loop effective action T is given by the expression 

iT = ~ Tr In K maP + Tr In Mj - ^ Tr In C afi (5.5) 

Notice that we use a different definition of M and that is why the sign of Tr In C in (B.5) is 
opposite to that of expression (|T8|). If we denote by X the derivative of the quantity X with 
respect to a generic parameter introduced by the gauge fixing functional, then the derivative 
of (B.5) with respect of such a parameter can be estimated with the use of (B.4). The result 
is 

^X-y ^ja-y _|_ ^ (J ^X^ 1 

The difference between the one-loop corrections in two different gauges, T(a) — T(a) is ob- 



•t = ~\ Tr G„ ^- ^ 

2 dg p(7 bgxr 



(5.6) 



tained by integrating (B.6) for all the gauge fixing parameters from a to a pB]| . The equation 



(B.6) gives the general form of the gauge parameters dependence of the one-loop effective 
action which is proportional to the motion equations, 5S/8g a p. For the higher derivative 
theory with iV > this implies the independence of the one loop divergent contributions on 
the gauge parameters a n ,7 n and (3 as proved in section 3. 

Appendix C 

The divergent part of the diagrams of Figure 3 can be calculated by means of the back- 
ground field method and Schwinger-DeWitt technique, using the generalization introduced 
in Ref. |57j and the methods developed in Ref. J29| (see also fl30"| ). In the present theory the 



calculation is technically simpler than the one carried out in p9[ . First, we notice that, by 
power counting, we only need to calculate the divergent terms of the form Z 2 where <fi is 
the background scalar, and Z is a divergent coefficient with dimension of mass 2 . Besides the 
splitting of the metric g pu — > g pu + h^ u into background and quantum fluctuations we have 
a similar splitting for the scalar field <ft — > <P + V 9 - We can therefore consider the background 
field as a constant. The second diagram contains purely gravitational loops, and it can be 
evaluated with the algorithm used in (|32[). Then one gets the result by replacing 



in the expression 



1^ (C.l) 
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The calculation of the divergences generated by the first diagram is more involved. Power 
counting tells us that this contribution only appears when all derivatives of the vertices of 
those digrams act on the internal lines. In the framework of the background field method 
this means that we need to keep only the terms with two derivatives of the quantum fields 
when performing the background field expansion of the £ -R 2 term in (|40| ) . On the other 
hand one can neglect all the non-leading order terms in the pure gravitational h — h sector 
and consider only the contribution of the higher derivative terms, because all the others 
do not contribute to the divergent counterterms. The same feature occurs for non-leading 
derivative terms in the scalar — sector. However the higher order terms in the mixed 
— h and h — <p sectors are relevant. In summary, it can be shown that all the divergences 
are contained in the expression (i/2) Tr In H where 



H 



and A,P,Q are given by 



P 



At 



fJtU 



v A v/ 



Q 



At 



a/3 



v A v r 



-2-° 



(C.2) 



.4 



Jfiu,af3 



P 



At 



£0 (5 A V 



a/3 



£0 {S > 



a/3 



9 XT 9a(3 



(C.3) 



This formula can be compared with the similar expression associated to fourth derivative 
gravity coupled to matter |30| . 



Thus, 



Tr In H = Tr In 



.4 



pu,a/3 





+ Tr In 







Q 



At 



a/3 





-1/2 


V A V T n- 3 



+ Tr In 



I 
V A V r D 




(CA) 



The first term in the last expression is the logarithm of the determinant of a c-matrix, 
and thus it is finite. The second term depends only on the metric but not on the background 
scalar, and therefore it can give only contributions to the renormalization of the gravitational 
sector. Since this renormalization has been already considered in the diagrams of Figure 2, 
we have to ignore here its contributions. Expanding the last logarithm in power series, one 
can keep only the terms with proper background dimension. Power counting tells us that 
only these terms are divergent (one can also use universal trace formulae of |37|] to check the 
convergence of the rest of the series) . The Z 2 -type divergence is given by the second term 
of the series and has the form 

1 r ~ A ■ F "At1 a (2) 

la/3 jXrpcr 



12 e 



d Xa 



-g 



p 



A 



fJ,V 



fiu,a(3 [Q^ ' 



(C.5) 
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Substituting the above expressions for A, P, Q and taking into account the contribution of 
the second diagram of Figure 3, we get the following divergent part of the effective action 
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FIGURE CAPTIONS 



Figure 1: Diagrams contributing to the renormalization of matter and gauge field (thick 
lines) interactions in a gravitational classical background (wavy lines). 

Figure 2: Diagrams contributing to the renormalization of gravitational, matter and gauge 
field interactions including quantum gravity corrections. 

Figure 3: Divergent diagrams contributing to the mass of scalar particles in the N — 1 
theory. The divergent contribution is generated by the non-minimal couplings 
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